Rescheduling trains by crossover tracks to promote service quality of urban rail transit under partial blockages

This study focuses on the rescheduling problem with disruptions that cause partial blockages in the urban rail transit (URT), contributing to extending the relative train rescheduling studies. The alternative driving measure (ADM), which could be regarded as one train rerouting measure, is used to skip the blocked section, and a mixed-integer nonlinear programming (MINLP) model is built based on it. Time-varying passenger flow as well as the turnaround process of rolling stocks is taken into consideration. To solve the model, a customized genetic algorithm is used to quickly generate high-quality solutions. Real-world data is studied and sensitivity analyses are taken to verify the feasibility and advantage of ADM. The results validate the proposed model and algorithm, as well as confirm that ADM shows significantly better performance than the practical operation measure in promoting passenger service quality of URT under partial blockage.


Introduction
Urban rail transit (URT) has been developing at a great high speed these years, which is broadly accepted for its large capacity, high service frequency, low energy consumption, and low pollution [1,2].As a result of its high passenger density and short headway, any unexpected technical or human incident may cause heavy disturbances to passengers' trips, pushing operators to focus on the passenger service [3,4].Cacchiani et al. [5] defined incidents as disturbances and disruptions in an overview paper, and rescheduling problems for them have been widely studied for decades and many achievements have been got.
Usually, there are several train rescheduling measures for the operator to choose from considering different situations.It is worth noting that though there are significant differences between URT systems and passenger railway systems, many rescheduling measures are similar for both systems, and researchers would use these measures for different systems by giving specific consideration to line conditions, technical standards, etc.Firstly, using the buffer time, which is also named retiming, could be regarded as an efficient way for real-time operation [6], which reschedules trains by adjusting the running time, dwell time, as well as the headway [7][8][9].Besides retiming, flexible stopping, which is also called skip-stop, is a well-known rescheduling measure that can speed up the circulation of trains and reduce stranded passengers [10,11].Both mentioned measures mainly focus on incidents that do not cause any blocked section in the line, while some disruptions could result in partial or complete blockages [12,13].As a substitute measure, bus bridging is gradually contended to be one important part of disruption management, which refers to using buses to transport passengers to their planned destination during metro disruption [14,15].Bus bridging shows a positive influence on relieving stranded passengers caused by disruption, but the passenger service quality inside the urban railway system could not be improved by this measure.
In previous studies which focus on rescheduling during a complete or partial blockage, the rerouting measure has been popular, which could be divided into passenger rerouting and train rerouting.Passenger rerouting means that passengers transfer to other trains or choose other routes to continue their journey to ensure they could reach their destinations during the disruption [16,17], which is mainly used in railway rescheduling or metro network rescheduling problems.Train rerouting could be used for both a network and a line [18,19], and for a line it contains two kinds of measures.One of them is short-turning, which has been widely studied in both railway and URT.Ghaemi et al. [20] considered short-turning trains in one direction at multiple stations adjacent to the blocked section, and Ghaemi et al. [21] investigated the impact of the disruption length based on the rescheduling measure which was shortturning.In URT, most of the rescheduling studies with the rerouting measure have used shortturning.Bes ˇinović et al. [1] proposed an integrated traffic management model for the complete blockages in urban railway lines, and Long et al. [16] studied the passenger-oriented rescheduling of trains and rolling stocks under disruption for a metro network.Li et al. [22] considered the short-turning to be one possible choice of the train turnaround and focused on the the resilience of urban rail systems.Short-turning could ensure the balance of the rolling stock and help improve service quality in cases of complete blockages.However, considering the partial blockage, passengers who want to travel for a long-distance crossing the blocked section can only wait for recovery with short-turning, which will significantly increase the travel time of these passengers, resulting in a decrease in passenger satisfaction.Worse yet, if these passengers are stranded for a relatively long time, they would prefer to choose other transportation modes to get to their destinations, which causes a significant decline in the service quality of URT.
The alternative driving measure (ADM), which is the other measure of train rerouting, shows significant strength in dealing with partial blockages by offering an efficient way to let passengers continue their trips during the disruption [23,24].Though very few, there are some scholars who have paid attention to taking ADM into account in the train rescheduling problem.Xu et al. [25] proposed a model with ADM during the partial blockage and verified the advantage of this measure in declining the delay of train services and ensuring the service balance at two terminals.However, their research didn't take passengers into account, while ADM could bring significant improvement in passenger service quality.Huang et al. [26] developed rescheduling models based on two strategies considering the partial blockage of the circular line, focusing on alleviating the inconvenience for passengers and regaining the nominal train regularity.One of the two used strategies in this research was the same as ADM.However, they didn't consider the rolling stock circulation, which is important under partial blockages, because it is difficult for rolling stocks in the blocked direction to finish their task and begin the new in time.Also, in their future research part, they mentioned one may extend those models to use the actual numbers of boarding and alighting passengers at each station, as they had done some simplification.
It is widely accepted that the train rescheduling problem with a large number of integer variables is challenging to get an optimal solution in acceptable solving time.Therefore, many researches applied heuristic algorithms as well as metaheuristics to get near-optimal solutions for the train rescheduling problems.Zhu et al. [27] studied the collaborative optimization of rescue operation and timetable rescheduling problem and solved the proposed model by a tabu search algorithm.Simulated annealing algorithm was used in Tamannaei et al. [28] for tackling the large-scale train rescheduling problems.Xu et al. [29] designed an efficient genetic algorithm to solve a passenger-oriented model for train rescheduling, and Han et al. [30] designed an algorithm based on the Non-dominated Sorting Genetic Algorithms-II (NSGA-II) to find the trade-off between two objectives of the proposed model.Zhang et al. [31] designed an iterative meta-heuristic for the general metro rescheduling algorithm to classify an accident and determine rescheduling solutions.For metaheuristics to solve the train rescheduling problem, some studies combined the rescheduled timetable with some other operation plans to be designed as the initial solution [29,30], while some other studies may have the operation plan to be the initial solution and use solvers such as Gurobi to get the rescheduled timetable [27].Then the iteration process would be taken to generate relative and better solutions, and finally find one near-optimal solution for the problem with some terminal conditions.
Table 1 shows the comparison of different publications of related rescheduling models and methods.For the names of columns, the time-varying passenger flow with OD is indicated by TPOD, and the turnaround process is represented by TP.
In summary, there are seldom studies using ADM as the rescheduling strategy for partial blockage in URT.In this paper, we have studied normal urban railway lines, i.e., acyclic lines, which could be regarded as the subset of a network if only one disruption occurs at one time, and ADM has been taken into consideration to continue all the services during the disruption.Compared with the few existing papers that apply ADM, this paper has contributions in considering the time-varying passenger flow with OD, as well as the turnaround process of rolling stocks.Specifically, the contributions and differences of this paper are given as follows: Blockage: no blockage is considered in the publication (N); partial blockage is considered (P); complete blockage is considered (C).

TPOD & TP: consider (C); non-consider (N).
Solution Method: model predictive control (MPC); the genetic algorithm (GA); the iterative algorithm (IA); the tabu algorithm (TS); the simulated annealing algorithm (SA); the integrated disruption management (IDM); the rolling horizon (RH); the capacity check algorithm (CCA). https://doi.org/10.1371/journal.pone.0296018.t001 1.This paper considers dynamic passenger flows, which are dealt with as actual numbers with origins, destinations, as well as certain arrival times, instead of arrival rates and alighting ratios for simplicity.Also, besides waiting and boarding, passenger behaviors are extended in this paper by introducing leaving after waiting for a relatively long time.Both improvements in the passenger flow are not considered in previous research which studied ADM in URT.
2. This paper develops a mixed-integer nonlinear programming model based on the ADM, which takes the turnaround process into account.Though the turnaround process contains an important part of rolling stock circulation, previous research using ADM in URT almost ignores it to simplify the rescheduling problem.A customized genetic algorithm is designed to meet the computing requirement of rescheduling, of which the population generation could produce effective child individuals, and a check process is used to ensure individuals' feasibility, which could help handle integer programming models with many complicated constraints by genetic algorithm.
The rest of this paper is structured as follows.In the Mathematical model section we describe the rescheduling problem and then propose a mixed-integer nonlinear model based on ADM.A customized genetic algorithm is designed in the Algorithm description section to solve the model, and in the Case study section we deal with the real-world data to validate the proposed model and algorithm, also verify the advantage of ADM in improving passenger service quality.At last, we deliver a conclusion of this paper and provide some further ideas in the Conclusions section.

Problem description
Basic descriptions.As aforementioned, incidents causing partial blockage would divide a bidirectional urban railway line into two double-line sections and one single-line section by two crossovers.Take the upstream direction to be the blockage direction as an example, and Fig 1 shows the three sections named as S 1 , S 2 , and S 3 respectively.Let the set of upstreamdirection stations S u J ¼ f1; 2; . . .; Ng, and the set of downstream-direction stations S d J ¼ fN þ 1; N þ 2; . . .; 2Ng.During the disruption period, the three sections could be regarded as five subsets of stations according to the two crossovers on both sides of the blocked section, and S u 1 ¼ f1; 2; . . .; Pg, S d . . .; 2N À Qg.Station P represents the last station before the blocked section in the upstream direction, and station Q is the station at the end of the blocked section.The shape of crossovers shown in Fig 1 is opposite, but other shapes are also feasible for this paper, only affecting some parameters for trains to pass.
To the best of our knowledge, for the commonly used measure in real-life operations for partial blockages in China, impacted train services in the upstream direction are mostly asked to wait at stations before the blocked section during the disruption, leading to unbalance of rolling stock circulation as well as passenger unpleasantness [25].After recovery, train services would be rescheduled by compressing the buffer time to decline delays, which is also called retiming, and finally achieve the original timetable.This commonly used measure is called practical operation, and stranded train is used in this paper to indicate the rolling stock stopping during the disruption.When the disruption occurs, train services in the upstream direction which have already passed station Q would not be affected, and it is similar for train services in the downstream direction which have already passed station 2N−Q.But for train services in the section between station P+1 and station Q, they would be stranded during the disruption.For these stranded trains, they are assumed to stop at the nearest station and wait for recovery.An example is given in Fig 2 to compare the rescheduling process of the two rescheduling measures, the practical operation measure and the ADM, in which the first train service in the upstream direction represented by the blue line is stranded at Station 3 due to the partial blockage.In the planned pattern of rolling stock circulation, the third train service in the downstream direction is connected with the first upstream train service.The blocked section in shows that during the disruption, the second upstream train service could continue its journey from Station 2 to Station 5 by the crossovers and downstream tracks, which is shown in the red line.Passengers in the first upstream train who want to travel to Station 4 or 5 could alight at Station 3 and board the second upstream train at the opposite platform, instead of waiting in the stranded train until recovery.Also, the third train service in the downstream direction is connected with the second upstream train service now, and it means the third downstream train service could begin its journey earlier with ADM than that with the practical operation method.In this way, ADM could show a significant advantage in maintaining rolling stock circulation and ensuring passenger reachability during the disruption.However, because the region from Station 3 to Station 4 is shared by two-direction train services as a single-line section, the capacity over this region would be reduced with ADM [26].Indeed, ADM is not suitable when the single-track section is too long, therefore, the blockage studied in this paper is considered to exist between two adjacent crossovers instead of a long range.
Descriptions of train services.In this paper, train service is defined as a rolling stock operating between two terminal stations.f is used as the index of train services in the upstream direction, and g is used for train services in the downstream direction.F and G represent the set of train services in each direction, respectively.j is used as the index of stations, j 2 S u J [ S d J , then j−1 and j+1 respectively indicate stations before and after station j along the train service direction.For train services in the upstream direction, there are two groups divided by whether the service passes the single-line section.Let F d be the set of train services in the upstream direction that needs to pass the single-line section, Namely, f 1 indicates the first train service in the upstream direction that needs to pass the single-line section and f k is the last one.f d is used as the index of train services that need to pass the singleline section, f d 2F d .
As aforementioned, the blockage doesn't exist for a long range, meaning that there would not be too many train services stranded in the blocked section.Therefore, this paper takes the situation in which there is only one stranded train service f 0 in the upstream direction stopping in the blocked section as an example.To find which service is the f 0 , we have where a o f 0 ;Pþ1 is the arrival time of service f 0 at station P+1 in the original timetable, and a o f 0 ;Q ; a o f 0 À 1;Q are similar; d o f 1 ;P is the departure time of service f 1 at station P in the original timetable; T 0 is the start time point of the blockage, and it is easy to find f 1 = f 0 +1, namely, f 1 is the first service practically behind f 0 .Formula (1) represents that f 0 has already entered the blocked section when the disruption occurs, while it has not left this section according to Formula (2).Formula (3) & (4) indicate that there is only one train service stranded due to the disruptions.j 0 is used to represent the last station served by f 0 , and for j 0 2 fP þ 1; P þ 2; . . .; Qg; a o f 0 ;j 0 � T 0 � a o f 0 ;j 0 þ1 .Correspondingly, we use g 0 to represent the last train service in the downstream direction to leave station 2N−Q before the disruption, and then g 1 indicates the first train service in the downstream direction that needs to pass the single-line section.For train services in the downstream direction we have We use G s = {g 1 , g 2 ,. ..,g h } to be the set of affected services in the downstream direction which need to pass the single-line section during the disruption, and g s is used to represent the index of these train services, g s 2G s .
Descriptions of rolling stock circulation.Considering the circulation of rolling stocks, it is important to ensure whether there is any backup rolling stock in the depot, and this would have a significant influence on the rescheduling process.Obviously, more backup rolling stocks could make the rescheduling problem under the partial blockage much easier.While in this paper, we take the limited amount of rolling stocks into consideration and only focus on the existing rolling stocks in the planned pattern, which could be a difficult but possible situation.There has already been a planned pattern for the connections of rolling stocks, which are described as given binary variables d o f ;g and d o g;f .There is d o f ;g ¼ 1 if service g in the downstream direction uses the same rolling stock of service f in the upstream direction, and d o g;f ¼ 1 means service f in the upstream direction uses the same rolling stock of service g in the downstream direction.Correspondingly, binary variables δ f,g and δ g,f are used to describe the connections of rolling stocks in the rescheduled pattern.
As the disruption occurs in the upstream direction, during the whole rescheduling period there is always Whether there is any difference between d o f ;g and δ f,g depends on the number of stranded trains caused by the blockage, and in the studied case with one stranded train, during the disruption there is and g2G.Then after the recovery, f 0 would continue its service at remained stations, d f 0 ;g ¼ d o f k ;g , and For the situation in which there is no train service stranded in the blocked section, it could be regarded as one simplified situation of the studied case.For all train services f2F, g2G, d g;f ¼ d o g;f and d f ;g ¼ d o f ;g .While for situations with more than one stranded train service, the rescheduling pattern for the connection of rolling stocks could be projected with the corresponding index, for example, d f 2 ;g ¼ d o f 0 ;g for two stranded trains, etc.All the situations mentioned have a similar rescheduling process except the pattern for the rolling stock circulation, and then we take the situation with one stranded train service as the case to be studied.
Descriptions of passengers.Procedures of passenger behaviors considered in this paper include arriving, waiting or leaving, boarding, and alighting, which is shown in Fig 3 .Parameter p jj 0 t represents the number of passengers that arrive at station j at time point t with the destination station j 0 , and p jj 0 t ¼ 0 when station j and j 0 could not form an OD pair.However, during the disruption, some OD pairs would be affected due to the blockage, in which there is at least one station, origin or destination, in the blocked section.For passengers who arrive at or expect to travel to station j in the blocked section during the disruption, j in their OD pairs should be replaced by 2N+1−j, which means the opposite platform of the same station.
Before the service arrives, there could probably be some leaving passengers because the waiting time is too long for them.Integer variable q j i is used to describe the number of passengers leaving station j for other transportation modes before service i arrives at this station, and P j i indicates the number of passengers finally board service i at station j.A parameter t O is given as a constant which depends on the results of a simple survey for the tolerance time for passengers to wait for one service in URT, which is related to q j i .Moreover, binary variable z j it ¼ 1 indicates that passengers arriving at station j at the time point t finally could board service i.In this paper, we take the off-peak period and the leaving behavior into consideration, and the limitation of vehicle capacity is ignored, meaning that whether passengers could board one service only depends on the time factor.

Problem assumptions
For rigour and simplicity, several assumptions are made for the problem as follows: 1.The number of used rolling stocks is a constant, namely, the disruption does not affect the technical condition of all rolling stocks and no backup rolling stock would be considered in the rescheduling process.Tracks and stations are all allowable for train services from both directions, and there is no difference among rolling stocks.
2. The stranded train would stop at the last served station during the disruption.Passengers on the stranded train would alight and go to the opposite platform for another service.Also, passengers who expected to arrive at the blocked station during the disruption would practically wait at the opposite platform to board one train service.
3. Disruptions occur in off-peak hours, and all passengers who have not left could board the first train service after their arrival without waiting for the next one.Passengers would choose to leave the URT for another transportation mode when the waiting time is longer than their tolerance time.

Model formulations
Based on problem description and problem assumptions, model formulations are given as notations, timetable and passenger flow formulations, as well as objective functions.Notations.All the relevant notations and parameters used in formulating the model are shown in Table 2.
Timetable formulations.1. Arrival and departure formulations In the rescheduled timetable, let a g,j and d g,j be the actual arriving and departure time at station j of service g in the downstream direction, there is: where w g,j is the actual dwell time of service g at station j, and r g,j−1 is the actual running time of service g from station j−1 to station j.Considering the practical operation requirement, there are limitations on the dwell time and running time: where w min j is the minimum dwell time at station j; r min jÀ 1 and r max jÀ 1 represent the limitations of running time between station j−1 and station j.In fact, the actual dwell time and running time would depend on the operation instructions and the limitations.
For train services in the upstream direction, there is: where w f,j is the actual dwell time of service f at station j; r f,j−1 is the actual running time of service f from station j−1 to station j; t c is the additional time of trains to pass the crossover and Eqs ( 13)-( 16) describe the process for train services f d entering and leaving the single-line section by the crossover.Especially, considering the stranded train f 0 and the last train service passing the single-line section f k , we have where T 1 is the end time point of blockage.Set of stations in the downstream direction Subsets of stations during the disruption according to It is particular for the arrival times of train services at the first station of each direction, as there are two possible cases.Take services in the upstream direction as an example, as services in the downstream direction are similar.Some of services in the upstream direction use rolling stocks which are not from services in set G, as their rolling stocks correspond to earlier services in the downstream direction which are not considered in the rescheduling process.Their arrival times at station 1 could be the same as the original timetable.While others which are connected with services g2G, the calculation should consider both the turnaround process and the departure time of service g at station 2N.Therefore, the arrival time of train service at the first station of each direction could be computed by where t b g is the time for the rolling stock of train service g to finish the turnaround task, t b f is the time for the rolling stock of train service f to finish the turnaround task, and both of them should satisfy the minimum turnaround time for rolling stocks, i.e.,

Headway formulations
To ensure the security of operation, consecutive train services should follow several headway limitations.Then we formulate the headway constraints for train service g in the downstream direction as where h 1 is the minimum headway for two consecutive train services to depart from one station; h 2 is the minimum headway for the following service to arrive at a station that the leading one departed from; d g−1,j is used to represent the time point of the last train service departed from station j before service g, and when g = 1, d g−1,j could be gotten from the original timetable.
Then for train services f in the upstream direction, there are where d f−1,j is used to represent the time point of the last train service departed from station j before service f.Also, when f = 1, d f−1,j could be gotten from the original timetable.

Alternative driving formulations
As aforementioned, train services in the upstream direction which need to pass the singletrack section during the disruption are indicated by f d 2F d , and F d = {f 1 , f 2 ,. ..,f k }.Also, g s 2G s correspondingly indicates services in the downstream direction.To ensure the security, there could not be both directions train services in the single-line section at the same time.Binary variables x g s ;f d and x f d ;g s represent orders for train services of entering the single-line section.x g s ;f d ¼ 1 means that after g s leaving the single-line section, it is f d in the opposite direction to be the next one entering the single-line section, and x f d ;g s ¼ 1 is similar.Then we have where t s represents the security headway for train services in different directions to enter the single-track section.After f k and g h , all train services run on the planned tracks due to the recovery of blockage.It is easy to find that in Eqs ( 30)-( 31), there are multiplication of variables, namely, these equations are nonlinear.Passenger flow formulations.There are several variables about the passenger flow as introduced in Problem Description, and they are all connected with the rescheduled timetable.Take service g as an example, and for service f the formulation is similar.t O indicates the tolerance time of passengers, then for passengers arriving between services g−1 and g, when a g;j À d gÀ 1;j > t O , there is z j g;t ¼ 1 if t � a g;j À t O ; if a g;j À d gÀ 1;j � t O ; z j g;t ¼ 1 is suitable for all passengers arriving at time point t for service g at station j.Eqs ( 33 ( where z j g;t is the binary variable to show whether passengers arriving at station j at the time point t finally could board service g.With z j g;t , we could formulate the number of passengers boarding service g and those leaving before service g as follows As aforementioned, when g = 1, d g−1,j could be gotten from the original timetable.Then for the total number of leaving passengers, we have Considering the leaving behavior of passengers, we use a variable ω to indicate the ratio of failed travel for passengers, namely, they could not board one service and have to leave the URT system.And for this ratio there is Objective functions.In this paper, passenger service quality is the total objective, and passenger cost is used to describe it.It is clear that the minimum passenger cost corresponds to the best passenger service quality, and the passenger cost could be considered by two parts due to whether passengers could board one service.For those who leave the URT system due to a relatively long waiting time, a parameter t p is given as the penalty time to calculate their costs.
Considering passengers who finally board one service instead of leaving, the time cost including average waiting and travel time is taken into account first.The average waiting time for passengers is widely accepted to analyze the service quality of URT, and the average travel time also has an important influence on passengers' traveling experience because people usually hope to arrive at their destination as soon as possible.For the sum of waiting time and travel time of all passengers there are: Eqs (39)-( 40) are nonlinear with variables to be the bounds of summary, and there are also multiplications of variables in these two equations.Then for the average waiting time and travel time of each passenger, we have Besides the traditional time cost, there is special cost of service deviation.The deviation of the rescheduling service pattern from the original one was mainly regarded as the cost of the company in previous research because it means more operation consumption.However, with the development of the Passenger Information System (PIS), passengers can get enough information in advance and have expectations about the time for services to arrive at or depart from one station.And when their expectations are not realized, they may feel uncomfortable.That is why this paper put the deviation of service time as one influential factor in passenger benefit.Since early arrivals and departures could be allowed in the URT, the total service deviation of the rescheduling pattern from the original one is given by With the absolute value operation, Eq (43) is also nonlinear.For each passenger, he or she would only board at most one service at one time without transfer, and we have the average deviation as where c represents the passenger cost, ω indicates the possibility for passengers to leave the URT, and the value of t p could be gotten from the practical operation of one specific line based on experience.Due to the nonlinear characteristics of objective functions as well as Eqs ( 30)-( 31), the proposed model ( 45) is a mixed-integer nonlinear programming model which is difficult to be solved by existing commercial solvers, and we have designed a customized GA to get the near-optimal solution in acceptable computing time.

Algorithm description
To meet the requirement of real-time rescheduling, efficient algorithms are required, such as MPC from the control-oriented, heuristic algorithms like the genetic algorithm (GA), etc. Heuristic algorithms are widely used due to the large scale of rescheduling problems, and a customized GA is designed in this paper to solve the proposed model, which could be more adaptable for this train rescheduling problem.Compared with the standard one, the customized GA could get an optimized solution in a short time.

Population generation
For each individual, the chromosome is composed of two parts, which are coded by service times in the two directions separately, as shown in Fig 4 .In this figure, cubes in blue represent times of services in the upstream direction without passing the single-line section, while yellow cubes represent those passing the single-line section in the upstream directions, as there are differences in the indexes of stations.Green cubes represent times of services in the downstream directions.All the cubes could be regarded as genes on each chromosome.
In the initialization process, solutions, which are shown as timetables of the rescheduling period, are generated by operating train services in chronological order under the limitation of traffic constraints shown in Section 'Timetable formulations'.In this way, the solutions of the initial group would be all feasible.Differences among feasible solutions are mainly formed by the flexible running time of each section and dwell time at each station, as well as the different turnaround time of each service.In the standard GA, the initialization process has more randomness in each individual which represents the rescheduled timetable, because train services could stop and wait as long as the operator requests.However, the large buffer time could cause large delay for the later train services of the solutions, which would result in worse solution quality.The customized algorithm sets upper bounds on dwell times and turnaround times to ensure the quality of these feasible solutions, and these bounds are calculated by adding some buffer time for the corresponding time in the original timetable.starts at 10:15 am and ends at 10:35, with the duration 20 min for the base scenario.The planned headway between train services of the original timetable is 240s, and the number of involved train services in each direction is 15.The number of passengers comes from the practical AFC data during the studied period.Other parameters of this base scenario are provided in Table 3, and the maximum and minimum running times of each section are given in Table 4.For the maximum running times, they are calculated by adding 10s to the planned running times of the original timetable, which is according to Wang et al. [33].While for the minimum running times, they are calculated by the section lengths and the maximum velocity of trains.Moreover, the minimum dwell time at each station is set as 20s.Specific data listed in Table 5 are the average value of ten experiments, which show the superiority of ADM in terms of a reduction of the passenger cost, namely, the improvement of passenger service quality, outperforming the practical operation measure by more than 13.5%.With applying ADM, there is a significant reduction in the service deviation, and also the ratio for passengers leaving URT declines.Each passenger would wait for nearly 10 seconds more than applying the practical operation measure, while his/her travel time is nearly 40 seconds lower.
Besides the function values, the number of leaving passengers at each station of the line is shown in Fig 10 .It could be found that applying ADM could significantly reduce the number of leaving passengers in the upstream direction, with the increased number of leaving  passengers in stations 16-23.But from the whole line viewpoint, the total number of leaving passengers with ADM is 3772, which is 28.80% lower than the number 5298 with the practical operation measure.While for the customized GA, a comparison is done in Table 6 to verify its efficiency in this rescheduling problem.Besides the standard GA, two classic metaheuristics are used to solve the proposed model, and they are the simulated annealing (SA) algorithm and the tabu search (TS) algorithm.Parameters of these algorithms are gotten from several experiments.The passenger cost, which is the value of the objective function, could indicate the solution quality of the algorithm with different improvements, and the computing time shows the efficiency of the algorithm to solve the proposed model.It could be found from Table 6 that SA shows the best computing time among the four algorithms, but it ends at the local optimal solution more easily.TS shows better performance in the solution quality compared with the standard GA as well as SA.However, the customized GA shows significant advantages in computational efficiency and solution quality compared with other three metaheuristic algorithms used for this proposed problem.Moreover, Fig 11 shows the convergence of the customized GA.

Sensitivity analysis
The coefficient of passenger feeling sensitivity.Considering the coefficient of passenger feeling sensitivity μ in the formulation (40), an analysis has been done to show how it affects the ADM rescheduling process with other parameters and conditions to be the same, and the result is shown in Table 7.The larger μ means that passengers are more sensitive to the deviation of timetable, namely, they could feel more unpleasant with the difference between the rescheduled timetable and the original one.
From Table 7, with the increased value of μ, the travel time and the passenger cost increase, and there is a decline in the deviation of the rescheduling service pattern from the original one.It could be found that after μ = 0.05, the decline speed of the deviation would slow down with the value of μ increasing.This may be because though passengers would feel more unpleasant, fewer deviations could be eliminated except those that couldn't be reduced due to the rescheduling constraints.

Headway of two consecutive trains
In this study, the headway of the original timetable is 240s, which is during off-peak hours and set as a parameter.the original timetable, comparing ADM with the practical operation measure.The optimization rate means the relative difference between the passenger cost of ADM and that of the practical operation measure, i.e., it is calculated by "(the passenger cost of the practical operation measure-the passenger cost of ADM) / the passenger cost of the practical operation measure".It could be found from Fig 12 that for each analyzed headway, ADM shows more than 13% advantage in reducing the passenger cost compared with the practical operation measure, meaning that ADM is better than the practical operation method with these original headways.For the basic scenario in which the disruption duration is 20 min, the optimization rate of ADM grows with the increase of planned headway, showing that the proposed method may have more advantage on the original timetable with more buffer time.

Disruption duration
In this paper, disruption duration is taken as one given parameter.Then a comparison is made between the passenger cost of ADM and that of the practical operation measure with different lengths of the disruption period, and the results are shown in Table 8.From Table 8 and Fig 13, with the increased disruption duration, the passenger cost of both ADM and the practical operation measure would increase, as well as the increasing optimization rate, meaning that ADM shows better efficiency with a longer disruption duration.When the duration is less than 600s, namely, for disruptions under 10 minutes, there is little difference between ADM and the practical operation measure, as the optimization is less than 5%.

Tolerable waiting time
Tolerable waiting time, which should be a variable resulting from personal character and the travel purpose in the practical research, is simplified in this study as a constant.And with the scheduled headway of 240s, the tolerable waiting time is analyzed from 240s to 360s in Table 9, and the variation of passenger cost is shown in Fig 14.
From Table 9 and Fig 14, the increase in the tolerable waiting time leads to a significant decrease in the leaving ratio, and the passenger cost also declines with the increased tolerable waiting time, which means improvements in passenger service quality.Tolerable waiting time could be further studied by taking the passenger difference into account, and the analysis of  the tolerable waiting time could contribute to the work of leading passengers under disruptions.

Parameters for the algorithm
As mentioned in Table 3, there are several parameters set for the proposed algorithm.Experiments are taken to analyze different combinations of parameters for GA, and the chosen parameters are the population size R, the mutation probability p m as well as the given value for termination ε. Results of different combinations are shown in Table 10 as follows.
From Experiments 1-3, 4-7 and 8-11, it could be found that with other parameters to be the same, when the given value for termination ε is larger, the passenger cost would be larger and the computing time would be less.For the experiments of which ε = 0.01%, the computing time is acceptable with better solutions, while Experiments 1, 4, 8 which have ε = 0.005% could get relatively better solutions than other ε with quite longer computing times.Therefore, we prefer to using ε = 0.01% in the case study as well as other experiments.Also, results of these experiments show that the lower mutation probability would lead to local optimal solution more easily for the proposed problem.This may be because the applied crossover strategy would choose one train service to be exchanged, while for younger children it would be more possible for them to have same times of some train service, namely, crossover may make no changes at all.That may be why larger mutation probability could help to avoid getting local optimal solution early.For Experiments 12-17 with the same ε = 0.01%, different population sizes are analyzed with different mutation probability.The results show that larger population size would cause more computing time, while smaller population size would get local optimal solutions within significantly short time.Because larger population size means that each generation more feasible solutions need to be generated, which requires much more computing

Conclusions
In this paper, a mixed-integer nonlinear model is developed with the rescheduling measure ADM to solve the problem with a partial blockage caused by the disruption.The passenger cost is used as the objective function, which quantifies the passenger service quality of URT.The model is solved by a customized GA to get an optimized solution in a short time, and realworld data of one subway line in China is studied in the case study to verify and analyze the advantage of ADM in promoting passenger service quality compared with the practical operation method.The results validate the proposed model and algorithm, and it is confirmed by a sensitivity analysis that with different original headways during off-peak hours, ADM is always feasible for rescheduling under partial blockage and better than the practical operation method.For the determined disruption duration, ADM has the most effective performance at one specific original headway, and for headways longer than it, the optimization rate would decline for the better performance of the practical operation method.Other sensitivity analyses have been studied, and the result of analyzing the disruption duration shows that ADM may be more adaptive for the scenario with long disruption durations.The effect on improving the passenger service quality of ADM is positively related to the tolerable waiting time of passengers, but negatively related to the sensitivity of passenger feeling about the deviation.
The presented work contributes to extending the research using ADM to deal with partial blockage in URT, and compared with related studies, time-varying passenger flow and the turnaround process are firstly taken into consideration in this paper.Nevertheless, it could be studied to integrate short-turning and ADM as intelligent rescheduling measures for disruptions, and the rescheduling pattern for each train could be more flexible in the real-time rescheduling process.Furthermore, the tolerable waiting time could be set as a variable for different groups of passengers and combined with passenger behavior analysis.Moreover, this work only considers the application of ADM to normal urban railway lines.However, it could be a new challenge to adopt such kind of flexible routing to the train rescheduling problem in URT network, since the location of disruption and the impact of transfer passenger flows need to be incorporated into the current model.Also train services in different lines may share their rolling stocks in a network problem, which is another interesting research direction in the future.

Fig 1 FNmin jÀ 1 =r max jÀ 1
Set of train services in the upstream direction, |F| is the total number of train services in the upstream directionGSet of train services in the downstream direction, |G| is the total number of train services in the downstream direction F d Set of train services in the upstream direction that need to pass the single-line section during the disruption, F d �F G s Set of train services in the downstream direction that need to pass the single-line section during the disruption, G s �G Indices Definition j Index of stations, j 2 S u J [ S d J Index of the last station in the upstream direction P Index of the last station before the blocked section in the upstream direction Q Index of the last station of the blocked section in the upstream direction f Index of train services in the upstream direction, f2F g Index of train services in the downstream direction, g2G f d Index of train services in the upstream direction that need to pass the single-line section during the disruption, f d 2F d g s Index of train services in the downstream direction that need to pass the single-line section during the disruption, g s 2G s f 1 Index of the first train service in the upstream direction that needs to pass the single-line section f k Index of the last train service in the upstream direction that needs to pass the single-line section f 0 Index of the stranded train in the upstream direction g 1 Index of the first train service in the downstream direction that needs to pass the single-line section g h Index of the last train service in the downstream direction that needs to pass the single-line section g 0 Index of the last train service in the downstream direction to leave station 2N−Q before the disruption Parameters Definition d o f ;j =a o f ;j Departure or arrival time of service f at station j in the original timetable d o g;j =a o g;j Departure or arrival time of service g at station j in the original timetable w min j Minimum dwell time for train services at station j r Minimum/maximum running time for train services between station j−1 and station j d o g;f Given binary value, d o g;f ¼ 1 if service f in the upstream direction uses the same rolling stock of service g in the downstream direction d o f ;g Given binary value, d o f ;g ¼ 1 if service g in the downstream direction uses the same rolling stock of service f in the upstream direction δ g,f Binary value, δ g,f = 1 if service f in the upstream direction uses the same rolling stock of service g in the downstream direction δ f,g Binary value, δ f,g = 1 if service g in the downstream direction uses the same rolling stock of service f in the upstream direction t c Additional time for train services to pass the crossover t b min Minimum time for one rolling stock to finish the turnaround task p jj 0 t Number of passengers that arrive at station j at time point t with the destination station j 0 t O Tolerance time of passengers h 1 Minimum headway for two consecutive train services to arrive at or depart from one station (Continued )

Fig 7 .
Fig 7. The original timetable without the disruption.https://doi.org/10.1371/journal.pone.0296018.g007 which need to pass the single-line section during the disruption.In both Figs 8 and 9, the blue lines indicate upstream train services which don't need to pass the single-line section, and the green ones show train services in the downstream direction.From Fig 8, train services in opposite directions enter the single-line section alternately during the disruption, which offers a balanced service for passengers in both directions.Moreover, the dotted lines in the two figures represent planned services in the two directions.

Fig 10 .
Fig 10.The number of leaving passengers at each station.https://doi.org/10.1371/journal.pone.0296018.g010 Fig 12  shows the variation of optimization rate with different headways of

Table 2 .
(Continued) g /x g,f Binary variable, x f,g /x g,f = 1 if service g/f enters the single-line section when service f/g in the opposite direction leaves the section j g;t ¼ 1 if passengers arriving at station j at the time point t finally could board service f/g https://doi.org/10.1371/journal.pone.0296018.t002